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Abstract 

We consider the approach to self-similarity (or dynamical scaling) 
in Smoluchowski's equations of coagulation for the solvable kernels 
K(x, y) — 2, x + y and xy. In addition to the known self-similar 
solutions with exponential tails, there are one-parameter families of 
solutions with algebraic decay, whose form is related to heavy-tailed 
distributions well-known in probability theory. For K — 2 the size 
distribution is Mittag-Lefher, and for K = x + y and K = xy it is a 
power-law rescaling of a maximally skewed a-stable Levy distribution. 
We characterize completely the domains of attraction of all self-similar 
solutions under weak convergence of measures. Our results are analo- 
gous to the classical characterization of stable distributions in proba- 
bility theory. The proofs are simple, relying on the Laplace transform 
and a fundamental rigidity lemma for scaling limits. 
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1 Introduction 

Smoluchowski's coagulation equations provide a mean field description of 
several processes of mass aggregation in nature. We study the evolution 
of n(t,x), the number of clusters of mass x per unit volume at time t. 
Clusters of mass x and y coalesce by binary collisions at a rate governed by 
a symmetric kernel K(x,y), whence 



^r(^) = o / K(x-y,y)n(t,x-y)n(t,y)dy- K(x,y)n(t,x)n(t,y)dy. 



All microscopic interactions are subsumed into the agglomeration rate kernel 
K, and the process is assumed to be stationary in space. A broad survey of 
applications, especially in physical chemistry, may be found in the article by 
Drake [14]. Equation (1.1) has been used in an amazingly diverse range of 
applications, such as the formation of clouds and smog [17], the clustering 
of planets, stars and galaxies [32], the kinetics of polymerization [35], and 
even the schooling of fishes [29] and the formation of "marine snow" (see 
[21]). In the past few years, there has been a resurgence of mathematical 
interest in the field, largely due to the work of probabilists. An influential 
survey article by Aldous summarizes the recent state of affairs [4]. 

An issue of importance for homogeneous kernels, kernels that satisfy 
K(ax,ay) = aP K(x,y), is the phenomenon of dynamical scaling for all ini- 
tial data in a universality class. Mathematically, this corresponds to the 
problem of existence of scaling or self-similar solutions and characterization 
of their domains of attraction. In the aerosols community the relevant rubric 
is the theory of self-preserving spectra, which is treated at length in Fried- 
lander's book [17] and the extensive survey of Drake [14]. For a large class 
of kernels with 7 < 1, there is numerical evidence that solutions evolve to 
a self-similar form [23]. There are also physical self-consistency arguments 
that have been used to derive asymptotics for scaling solutions [34]. In the 
case 7 > 1, for a general class of kernels it is known that solutions must lose 
mass (presumably to infinite-mass clusters) in finite time [15, 19], but there 
is no general rigorous result on the precise nature of this blow-up in mass 
transport. In several instances the known solutions have unphysical diver- 
gences such as infinite mass. Thus, a general existence theory for finite-mass 
self-similar solutions, for example, would be of some value. 

The kernels K(x,y) = 2,x + y and xy play a special role, as (1.1) can 
then be solved by the Laplace transform. It is widely known that each 
of these kernels admits a self-similar solution with exponential decay (see 
Table 2 in [4]). These kernels are also special since certain solutions to (1.1) 





(1.1) 
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can be viewed as ergodic averages in beautiful probabilistic constructions, 
involving thinning of renewal processes {K = 2), and tree- valued Markov 
processes and their self-similar limits (K = x + y,xy) [2]. The additive 
kernel also figures in interesting recent applications given by Bertoin. It 
provides a natural probabilistic interpretation of a sticky particle model 
related to Zeldovich's model of gravitational clustering [8]. Also, the known 
self-similar solution appears in a simple model of turbulence, the inviscid 
Burgers equation with Brownian-motion initial data, as the characteristic 
measure for a Poisson point process which describes the shock strengths [6]. 

In short, aside from heuristics and numerics, there are no rigorous math- 
ematical proofs of the existence of self-similar solutions and the approach to 
self-similar form for general kernels (see [4, Sec 2]). And there are only a few 
partial results for the solvable kernels: For K = 2, Kreer and Penrose [22] 
proved local uniform convergence to the scaling solution under some techni- 
cal hypotheses on initial data. (Also see [11] regarding the discrete case.) A 
simple weak convergence theorem in this case follows from a classical result 
on the thinning of renewal processes [4]. In a recent article, Deaconu and 
Tanre proved a weak convergence result for all three kernels, but under re- 
strictive hypotheses on initial data [12]. Aldous and Pitman have studied the 
"eternal additive coalescent," and Bertoin has characterized "eternal solu- 
tions" to the Smoluchowski equation with additive kernel, solutions defined 
globally for — oo < t < oo [3, 7]. Bertoin showed these solutions correspond 
to the Levy measure of a first-passage process related to Levy processes with 
no positive jumps, and as a particular consequence he derived a new family 
of self-similar solutions related to the Levy stable laws of probability. 

In this article we find new families of self-similar solutions for the con- 
stant kernel, rederive the self-similar solutions to the additive and multi- 
plicative kernels by analytical means, and characterize all possible domains 
of attraction under weak convergence for all the solvable kernels. We show 
that 

1. For each of the solvable kernels, Smoluchowski's equation admits a 
one-parameter family of scaling solutions parametrized by a number 
p £ (0, 1] that characterizes the rate of divergence of the (7 + l)-st 
moment of the number density. For p = 1 these solutions reduce 
to the known solutions with exponential tails, while for < p < 1 
the number density has algebraic decay ("fat tails"). For K = 2 
(7 = 0) the normalized size distribution is a Mittag-Leffler distribution 
as studied by Pillai [31]. For K = x + y (7 = 1) and xy (7 = 2) the 
7-th moment distributions are transformed by power-law rescaling to 
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the Levy stable laws of probability theory (see 6.5 and [7]). 

2. The domains of attraction (under weak convergence of measures) for 
any scaling solution is determined by a condition on the tails of the 
initial data — the algebraic rate of divergence of the (7+ l)-st moment. 
A precise characterization is given via Karamata's notion of regular 
variation. In particular, with suitably diverging (7 + l)-st moment 
there are initial data for which there is no convergence to any self- 
similar solution. 

The self-similar solutions can all be captured by expressing their 7-th mo- 
ment distribution in the general form 

x 7 n(t, x) = m 7 (t)A 7 (t)- 1 / Pi7 (xA 7 (t)~ 1 ), (1.2) 

where explicitly, with (3 = p/(l + p), 

m {t)=t-\ mi(t) = l, m 2 (t) = (l-t)-\ (1.3) 
\ (t)=t 1 /o, \i(t)=eW, A 2 (t) = (l-t)- 1 /' 3 , (1.4) 

and the f pn are probability densities given by 

00 f_1 \k-l pk-l 

f P ,i(x) = f P ,2(x) = ~E J r(l + fc-fc/3)sinfe7r/3. (1.6) 

71 k=i 

We work with measure valued solutions to (1.1) denoted is t , where u t ((a, b]) 
denotes the number of clusters with size x G (a, b}: Vt{(a, b]) = J^n(t, x) dx 
if vi has integrable density n(t,x). For each of the solvable kernels, we 
associate a natural probability distribution function F(t, x) to the solution: 

rx I pod 

F(t,x) = ^ y^t(dy) J ^ y^ t (dy) ■ (1.7) 

This is the size-biased distribution for K = 2, the mass distribution for 
K = x + y and the second moment distribution for K = xy. We are inter- 
ested in necessary and sufficient conditions for the convergence of a rescaling 
F(t,X(t)x) to a nontrivial limit F*(x). Our results may be summarized in 
the following. 
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Metatheorem. For the kernels K(x, y) = 2, x + y, xy with degree of ho- 
mogeneity 7 = 0, 1,2 respectively, let T 7 = oo for 7 = 0, 1 and T 7 = T ge \ 
for 7 = 2. Then for any solution of Smoluchow ski's coagulation equation, 
there is a rescaling \(t) and a nontrivial probability distribution function 
such that 

lim F(t, X(t)x) = F*(x) at all points of continuity of F* 
if and only if 

X 

y 1+1 vo(dy) ~ x l ~ p L(x) as x — > 00, 

where p £ (0, 1] and L(x) is a function slowly varying at infinity. In the 
converse implication, F* must be a rescaling of F pri {x) = f£ f pa (y)dy. 

Precise statements are deferred to Theorem 5.1, Theorem 7.1, and The- 
orem 8.1. We write the results in this form to stress the analogy with 
the classical characterization of the Levy stable distributions in probabil- 
ity theory [16]. For the additive and multiplicative kernels the analogy is 
an intimate relation with distributions for asymmetric Levy flights: the self- 
similar solutions can be transformed by a power-law rescaling into maximally 
skewed a-stable Levy distributions. A deeper understanding of our results is 
obtained from Bertoin's study of eternal solutions [7]. The eternal solutions 
are analogous to infinitely divisible distributions of probability theory. Let 
vi be the value of an eternal solution at time t. Then, loosely speaking, for 
any s < t the measure v s decomposes v t such that u t is reconstituted from 
v s under coagulation. This heuristic statement is made precise by Bertoin's 
characterization of Levy pairs (a 2 , A) for the eternal solutions. This is the 
analogue of the classical Levy-Khintchine characterization of the infinitely 
divisible distributions [16]. Among the infinitely divisible distributions, the 
stable distributions are of special interest, and their Levy canonical mea- 
sures are pure power laws. And indeed, the self-similar solutions to Smolu- 
chowski's equation (1.1) have Levy pairs corresponding to pure power laws: 
a 2 = 0,A(dx) = cx~ a ~ l dx, 1 < a < 2 and a 2 = 1 and A = for a = 2 
exactly as in the classical characterization. When viewed in this context, 
our theorems are entirely natural. 

For K = 2 the main theorem may be interpreted probabilistically as 
a stability result for renewal processes on the line under uniform thinning 
(see [24]). For K = x + y the results are related to Burgers turbulence, for 
solutions of the inviscid Burgers' equation when the initial velocity is given 
by a Levy process with no positive jumps [6]. 
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We exploit the analogy with the classical limit theorems of probabil- 
ity to obtain simple proofs of optimal theorems. The proofs involve little 
more than the solution for the Laplace transform and a fundamental rigidity 
lemma that characterizes scaling limits via functions of regular variation [16, 
VIII. 8. 3]. And the analogy extends much further. The central limit theorem 
is perhaps the most intensively studied result in probability theory. Thus, 
we can demand stronger forms of convergence as in expansions related to 
the central limit theorem. In companion articles we plan to study (a) uni- 
form convergence of densities to the self-similar solutions with exponential 
tails (in analogy with the uniform convergence of densities in the central 
limit theorem) [27], (b) metric estimates (in analogy with the Berry-Esseen 
theorem) and (c) large deviation estimates. We have found proofs for (a), 
and partial results for (b) and (c), that follow easily from a combination 
of the solution formula and the classical method of characteristic functions 
outlined in Feller [16]. 

It is worth remarking that the folklore in the applied literature is that 
the scaling solutions are unique (e.g., see [34, 17]). This is false in gen- 
eral (though the solutions with exponential tails are indeed special — they 
attract all solutions with finite (7 + l)-st moment). With hindsight, this 
nonuniqueness is not surprising. The existence of a one-parameter fam- 
ily of scaling solutions is well known in physically related mean field models 
which show coarsening, such as the Lifshitz-Slyozov- Wagner model [28] , one- 
dimensional models for the coalescence of droplets [13], and in cut-and-paste 
models of coarsening [18]. From the mathematical point of view the funda- 
mental role of regular variation in branching processes is well established [9] 
and it is only natural that it should reappear in the "dual process" of coales- 
cence. We conjecture that analogous results hold for general homogeneous 
kernels, but these lie beyond our techniques based on the Laplace transform. 

2 Well-posedness for measure valued solutions 

2.1 Desingularized Laplace transforms 

Smoluchowski's equations determine a process of mass transport, and mea- 
sure valued solutions are an appropriate mathematical abstraction that con- 
tains solutions to the discrete and continuous coagulation equations within 
a unified framework. Norris recently proved several strong results for well- 
posedness of measure valued solutions, but these do not apply with quite 
the generality we prefer for K = x + y and K = xy [30]. We work with 
a somewhat different notion of solution motivated by the explicit solution 
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obtained with the Laplace transform. The use of the Laplace transform for 
these kernels is classical [14], and aside from trivial changes of notation, 
many of the equations below may be found in Bertoin's article [7]. While 
our primary aim in writing this article is not to tackle the question of well- 
posedness, we show below that the divergence of the self-similar solutions 
necessitates some care in the definition of solutions. The impatient reader 
may skim through this section making a note of the main theorems and 
the explicit solution formulas. Our primary source for background on the 
Laplace transform is Feller [16]. 

Let E denote the open interval (0, oo) and let A4+ denote the space of 
positive Radon measures on (0, oo). We interpret the number of clusters of 
size x 6 (a, b] per unit volume as v((a,b\) for u € M.+. We use the same 
letter to denote the distribution function of the measure, writing v(x) = 
v((0,x]), if this quantity is finite. We let m p = f E x p v(dx) denote the p-th 
moment of the measure, so uiq is the total number of clusters and mi is the 
total mass. 

We let r)(s) be the Laplace transform of v{x) defined by 



The last equation holds when v has a density n. In what follows we need to 
work with time-dependent measures v t for which the total number of clusters 
and/or total mass may be infinite. Consequently, it is more convenient to 
work with the variables (the "desingularized Laplace transform" ) given by 



The variable u = d s ip has the important physical interpretation that it is 
the Laplace transform of the mass measure. Probabilists will recognize the 
obvious similarity to the Levy-Khinchine representation. The equations of 
evolution in terms of these variables are extremely simple: 



V (s)= f e' sx u{dx)= f e' sx n{x)d. 



JE J E 



'X. 




(2.1) 



d t if = -if 2 
d t ip - fd s f = -Lp 
d t ip - ipd s ip = 



for K = x + y 
for K = xy. 



for K = 2 



(2.2) 
(2.3) 
(2.4) 



We will construct measures using these equations and establish that they 
are solutions of Smoluchowski's equation in an appropriate weak sense. 
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The function 1 — e~ sx does not have compact support in E, but has 
finite limits at and oo. A simple way to treat these limits is to consider 
E = [0, oo], and consider continuous functions on E, where /(oo) always 
means lim^oo f(x). 

Definition 2.1. 1. C(E) is the space of continuous maps / : E — > R 
equipped with the norm ||/||c<(b) = su Px 

2. C k (E) consists of k times continuously differentiable functions on E 
such that /,..., G C(E). It is equipped with the norm ||/||c<fc(^) = 

\\f\\c(E) + --- + \\f {k) \\c(Ey 

3. C k (E) is the subspace of C k (E) with compact support in E. 

4. £fc denotes the subspace of C k (E) of functions whose derivatives up to 
order k decay exponentially, that is X^j=i \f^\ x )\ — Cje~ ax for some 
a > 0, and whose derivatives up to order k — 1 vanish at zero, that is 
/O')(0) = for j < k. 

The following classical approximation lemma shows that the functions 
1 — e~ sx span a dense set in C(E). 

Lemma 2.2. (a) Let f G C(E). Then for every s > there is a sequence 
Pn(x; s) = Y2=i a n ,k{s){l - e ~ sx ) k such that 

lim \\P n (x;s)- f(x)\\ c{ E) = 0. 

(b) If f G £k then we also have Wran^^ \\P n (x; s) — f(x)\\ C k^ = for 
sufficiently small s. 

Proof. The problem may be reduced to polynomial approximation on the 
unit interval by the transformation y = 1 — e~ sx , and g(y) = f{x). Then 
g G C([0, 1]) for / G C{E). Assertion (a) now follows from Weierstrass' 
approximation theorem. A particularly useful choice are the Bernstein poly- 
nomials B nt g(y) of g, and P n (x;s) = B n , g {y). Suppose g G C k [0, 1]. It is 
then classical that lim^oo \\g(y) — -E>n,3(y)||c* fc [o,i] = [25, p. 25]. Thus, in 
order to obtain (b) it suffices to show that g G C k [0, 1]. A decay assumption 
on / is warranted, because by the chain rule 

sf(y) = f'(x)—, g"(y) = ^(f"(x) + s fXx)) 1 etc. 
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But / G £ k , so that £*=i |/^(a;)| < C/e"^. Hence, for s < a/k we have 
lim^i </(y) = . . . = lirrij^i g^ k \y) = and g £ C k [0, 1]. Thus, given any 
e > 0, there is an n(s) such that \g'{y) — B' ng (y)\ < e. The change of 
variables now works to our advantage, for we have \f'(x) — Pn(x)\ < ese~ sx . 
A similar calculation holds for all k derivatives proving (b). □ 

2.2 A weak formulation for measure valued solutions 

Following Norris, we will generalize Smoluchowski's equation as follows. To 
every finite, positive measure v we associate the measure L(y) defined by 
duality with continuous functions with compact support. 

(f,L{v))= l -l f [f(x + y)- f(x)- f{y)]K{x,y)u{dx)u{dy). (2.5) 

J E J E 

It is then natural to consider the weak formulation 

(/, u t ) = (/, vq) + f\f, L(v T )) dr, for every / G C C {E). (2.6) 
J o 

This suffices for the case K = 1, but it is insufficient for K = x + y and 
K = xy. The self-similar solutions to these kernels are not finite measures, 
and consequently they are not solutions in the sense of Norris, since they 
fail condition (3) in his definition [30, p. 80]. 

The basic obstruction is that L{y t ) is not a measure in general, since 
{f,L(ut)) may not be finite for all continuous functions. The reason is that 
even though / may have compact support in E the function Tf(x,y) := 
f(x + y) — f(x) — f(y) does not have compact support in E x E, and may 
not be integrable with respect to the product measure v ® v. 

Here is a counterexample for K = x + y. Let x( x ) De t ne indicator 
function for the interval (0, 1) and let v(dx) = x~ 3 / 2 x( x )dx + 5(x — 2). 
Let / be a continuous function with support in [a, b] = [2,5]. Then the 
values of Tf are as shown in Figure 2.2. Notice that v((x, oo)) diverges like 
0{x^ 1 / 2 ) as x —> 0, but v has finite mass. Thus, in order for the integrals 
in the definition of L to converge, it is necessary that there be suitably 
rapid cancellations as we approach the boundaries. One may show that the 
integral is finite on all regions except near the axes in the shaded regions. 
However, we explicitly compute that for any small 5 > 0, 

/ / [f(x + y)-f(x)-f(y)](x + y)v(dx)v(dy) 
Jo Jo 

5 

f{2 + y){2 + y)y-^ 2 dy. (2.7) 
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This is evidently infinite if f(x) rises sufficiently steeply for x > 2. Therefore 
L(v) is not a measure. 

This means that the space of continuous functions is not appropriate as 
a space of test functions in (2.6). The smaller spaces <? 7 serve as a suitable 
substitute. 

Definition 2.3. For each kernel K(x,y) = 1, x + y, xy with degree of 
homogeneity 7 = 0, 1, 2 respectively, let T 7 = oo for 7 = 0, 1 and T 7 = T ge \ 
for 7 = 2. We say that a map t ^ u t '■ [0, T 7 ) ^ M+ is a solution to 
Smoluchowski's coagulation equation if 

1. m 7 (0) = J E x 1 VQ{dx) < 00. 

2. For all compact sets B C E, the map 1 1-> vt{B) is measurable. 

3. Jq m 7 (r) 2 dr < 00 for all t G (0, T 7 ). 

4. For all / <E £ 7 and f G [0, T 7 ) we have 

(/,^) = (/^o)+ f\f,L(u T ))dr. (2.8) 

JO 
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2.3 Existence and uniqueness for the constant kernel 

We will set K = 2 instead of the usual convention of setting K = 1, since 
it simplifies several calculations (we actually revert to an older convention, 
see for example, eq. 459 in [10]). 

Theorem 2.4. Let vq G M.+ be a finite measure. Then Smoluchowski's 
coagulation equation with kernel K = 2 has a unique solution with initial 
data vo, and this solution is determined by the solution of (2.2). 

Theorem 2.4 is a consequence of [30, Thm 2.1]. We will prove it anew 
with the Laplace transform, as the explicit solution formula is needed later. 
Let ft denote the number distribution at time t, and (p(t, •) be determined 
from v = vt by (2.1) for each t. Then formally <p should solve the simple 
equation (2.2). For fixed s > 0, equation (2.2) is an ordinary differential 
equation with the solution 

p(t,s) = -^-. (2.9) 
1 + ipo{s)t 

Lemma 2.5. Assume vq G A4+ is finite. The formula (2.9) determines a 
weakly continuous map [0, oo) 3 1 1— > v t G M.+ with decreasing total number 
mo(t) = mo(0)/(l + mo(0)t). 

Proof. The solution ip(t, s) has the important property that its derivative 
is completely monotone for t > 0. This is because it may be written as a 
composition of positive functions with completely monotone derivative 

<XM) = ° (vo(s))- 

Recall that the derivative of (pit, s) is u(t, s) the Laplace transform of the 
mass measure, say fit(dx). Since u(t,s) is completely monotone, it follows 
that Utidx) G A4 + . Since u(t, s) is analytic in t, we see that the measures fj,t 
are weakly continuous by the classical duality between pointwise convergence 
of the Laplace transform and weak convergence of measures. That is for any 
continuous function / with compact support in E we have (/, fi T ) — > (/, fit) 
as t — > t. It follows that (/, v T ) — > (/, i/ t ) as r — > i where v t (dx) = x~ 1 fit(dx) 
is the number measure. The statement regarding mo(t) follows by taking 
s -> oo in (2.9). □ 

Remark 2.6. It is strange at first to consider the desingularized Laplace 
transform when the initial data is finite, and indeed the usual Laplace trans- 
form suffices. But (2.9) shows us that the solution is instantly regularizing 
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in the following sense. If ^o(°°) = ¥>o(°°) = °°, the solution satisfies 
lim <p(t,s) = lim — \ = -. 

s^oo s^oo 1 + t(po{S) t 

Thus, the number of clusters is finite for t > 0. Thus, (p(t, s) defines a 
natural solution even for an initially infinite measure. However, it is hard 
to verify (2.8) in this case (even for / 6 C^.(E)), and we restrict ourselves 
to finite measures in what follows. 

Proof of Theorem 2.4- Let us first show that the measures v t determined 
by the Lemma are a solution to Smoluchowski's equation in the sense of 
Definition 2.3. Property (1) has been assumed. It is easy to check property 
(2). The measures v% are weakly continuous. Thus for a fixed compact set 
B C E, the function 1 1— > ft(-B) is semicontinuous. Property (3) follows from 
Lemma 2.5. 

It is not a priori obvious that (1 — e~ sx , L(yt)) is indeed — ip(t,s) 2 . But 
the measures are finite since vt(E) < v$(E) < oo, and thus we may set 
/ = 1 — e~ sx in the definition of (/, L{y t )) to recover —ip(t, s) 2 . In particular 
this shows that (2.8) holds for / = 1 — e~ sx . This equation also holds if / is 
a monomial (1 — e~ sx ) k because 

(1 _ e~ sx ) k = ( k ) (-l) j e- jsx = -Y,( k -) t- 1 )^ 1 " e ~ jSX )' 

Given / G C(E) and e > 0, Lemma 2.2 guarantees an approximation P n 
with ||/ - P n || c( ^ < e. Then, \(f - P n , L(u t ))\ < 3e (u t (E)) 2 < 3e (ME)) 2 . 
Thus, 

(./>*>- <./>o> - [\f,L(u T ))dr 
Jo 

<(\f-P n \,u t ) + (\f-P n \,u )+ I \(\f-P n \,L(u T ))\dr 

Jo 

< e {ME) + ME) + 3t(v (.E)) 2 ) . 

This shows that the measures v t define a solution. 

Suppose vt and Ut are two solutions with the same initial data. Since 
/ = 1 — e~ sx G £q = C(E) and tp(t, s) = (/, u t ) < mo(t) for a.e. t, we can use 

(2.8) and part (3) of Definition 2.3 to obtain (2.3) in time-integrated form 
for each solution. It follows easily that for fixed s > 0, each ip(t, s) is C 1 in 
t and satisfies (2.3). But this equation has a unique solution <p(t,s) as in 

(2.9) . As we have noted in Lemma 2.5, (p(t,s) determines the measure v%. 
Thus vt = v t . □ 
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2.4 Existence and uniqueness for the additive kernel 

We always work with solutions of finite mass, but we do not assume that 
the number of clusters is finite. Therefore, L{y) will not be a measure. 
Nevertheless, it does define a bounded linear functional on the space of 
Lipschitz functions on E. 

Lemma 2.7. Let u G M+ with mi := J E xv(dx) < oo, and let L(y) be 
defined by (2.5) with K = x + y. Suppose f is Lipschitz and f(0) = 0. Then 

\(f,L(u))\<2mlUp(f). (2.10) 

Proof. By the symmetry of the integral in (2.5), we see that 

\(f,L(u))\< [ [ \f(x + y)- f(x)- f(y)\yu(dy)u(dx). 

J E JE 

The integrand is controlled by 

\f(x + y)- f{x) - f(y)\ < \f(x + y)- f(y)\ + |/(x)| < 2Lip(/)x. 
Thus we obtain 

\{f,L{v))\ < 2Lip(/) / / xyu(dy)u(dx) =2m?Lip(/). 

JE JE 

□ 

Theorem 2.8. Let vq G M.+ satisfy j E x v${dx) = mi < oo. Then Smolu- 
chowski's coagulation equation with kernel K = x + y has a unique solution 
with initial data vq, such that J E xvt{dx) = mi for all t G [0, oo). 

We will construct a solution using the desingularized Laplace transform 
and then prove its uniqueness. The evolution equation for ip is 

dt<p - fd s ip = -mnp. 

We may always normalize initial data such that m\ = 1, and we assume this 
in all that follows. Thus, we have 

dt<p - fd s ip = -Lp (2-11) 

It is striking that (2.11) is simply the inviscid Burgers' equation with linear 
damping. However, there is no shock formation, since the initial data are 
analytic with a completely monotone derivative satisfying d s (po(s) < 1. This 
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can be seen in the explicit solution below which is valid for all time. Since 
u = d s ip, differentiating (2.11) we have 

d t u — ipd s u = — u(l — u). (2-12) 

We solve (2.11) and (2.12) globally by the method of characteristics. Let 
s(t, a) denote the characteristic that originates at a at t = 0. Then we have 

ds dip du 

on a characteristic. We integrate (2.13) along the characteristics to obtain 

<p(t,s) = e-Vo(a), (2.14) 
s(t, a) - (p(t, s) = cr-ip (a), (2.15) 

U(t,S) = 77 7— T. (2.16) 

' 1 - (1 - e- t )u (a) v ' 

Lemma 2.9. Suppose vq{x) G Ai + with f E xv (dx) = 1. Then equation 
(2.11) determines a map [0, oo) 3 1 1— > vt & M+ such that 

1. f E xv t (dx) = 1 for all t. 

2. fit = xv t is weakly continuous. 
Proof. Observe that 

a - (po(a) = / (e~ sx - 1 + sx) i>o{dx) > 0, a > 0. 
Je 

Thus, by (2.14) and (2.15) 

s(t,a)=<T-(p (a)(l-e~ t )>0, a > 0. (2.17) 

The right hand side is a strictly increasing function of a for all t > 0, thus, 
the inverse map a(t, s) is well defined. Differentiating (2.17) with respect to 
s we find that 

da 1 , . _ t . . do" . . 

— = - --77 — t-t-tt, whence u(t, a) = e u (<t — . 2.18 

as 1 — (1 — e r )no(cr(s)) as 

Since uo is the Laplace transform of a positive measure, it is a completely 
monotone function of a. In order to show that u(t, s) is completely monotone 
in s, it suffices to show that da/ds is completely monotone in s (see Criterion 
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1 and 2 in [16, XIII.4]). We prove this as follows. Let us consider the 
sequence of iterates 

a Q (s) = s and <r n+ i(s) = s + (1 - e _t )<£> o (0"n(s)), n > 0. 

Clearly, |a n+2 (s) - a n+ i(s)\ < \a n+ i(s) - a n (s)\ so that o n {s) -> cr(s) the 
unique solution to (2.17). Moreover, we have 

= 1 + (1 - e £ )uo(o-n(s))- 



ds ds 

Thus itda n /ds is completely monotone, then so is da n+ \/ds. But dao/ds = 1 
is completely monotone. By induction, da n /ds is completely monotone for 
n > 1 and so is the limit da/ds. 

We may now conclude that the solution <p(t, s) to (2.11) defined by (2.14) 
exists for all t > 0, is unique, and has a completely monotone derivative 
u(t, s). u(t, s) defines a unique mass measure, say nt(dx). We see from the 
solution (2.16) that u(t, 0) = uq(0) = 1. Thus, the total mass f E fit(dx) = 
u(t, 0) = 1 for all t > 0. The measures ^ are weakly continuous since u(i, s) 
is analytic in time. □ 

Proof of Theorem 2. 8. Let us first check that the measures v% determined by 
(2.11) solve Smoluchowski's equation in the sense of Definition 2.3. Condi- 
tions (1) and (2) in Definition 2.3 are verified as in the proof of Theorem 2.4. 
Since v t has constant mass, it follows from Lemma 2.7 that the functionals 
L{y t ) are uniformly bounded on E\. In particular 

(l-e- 8X ,L{u t )) = -<p + <pd 9 tp (2.19) 

as desired. This shows that (2.8) holds for / = 1 — e~ sx , and thus for mono- 
mials (1 — e~ sx ) k . Given any / G £\ and e > we choose an approximation 
P n (x; s) as in Lemma 2.2(b) so that ||/ — P n \\c 1 (E) < £ - Then by Lemma 2.7 

\(f-P n ,L(u t ))\<2Up(f-P n )<2e. 

Similarly, 

\{f-P n ,v t )\ < ( \f(x)-P n (x;s)\u t (dx) <Lip(/-P n ) ( xv t {dx)<e. 

JE JE 

This shows that v t is a solution in the sense of Definition 2.3. 

Now suppose only that vt is a solution in the sense of Definition 2.3. The 
function / = 1 — e~ sx £ Si, and / < sx. For a.e. t > we have mi(t) < oo 
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and from tp = (f,v T ) < s _1 mi(t) it follows ip(t,s) is analytic in s with 
\dg(p\ < s 1 ~ k mi(t) for k = 1,2, . . .. Part (3) of Definition 2.3 ensures that 
for fixed s > 0, m\{t) is locally square-integrable on [0, oo) and (2.8) gives 
(2.11) in time-integrated form, since (2.19) holds for a.e. t and 

\(f,L(v t ))\ = \<pd 8 <p-<p\ < smi(*)(mi(*) + l). 

It follows that tp is continuous in t, uniformly for s in compact sets in E. 
Moreover, one can justify differentiating (2.8) in s and infer that d s (p is 
continuous. Then ip is a C 1 solution of (2.11) whence f(t,s) is uniquely 
determined by initial data. But as we have noted in Lemma 2.9, tp(t, s) 
uniquely determines the measure Vf Thus the solution is unique. □ 

2.5 Existence and uniqueness for the multiplicative kernel 

The multiplicative kernel differs from the constant and additive kernels, since 
it is not well posed for all time. But, the analysis can be formally reduced 
to the additive case by a change of variables. This is well-known [14], but 
we include it for completeness. 

The divergence of the classical self-similar solution is 0(x -5 / 2 ) as x — > 0. 
The total number and mass are infinite, but the second moment is finite. 
Therefore, we consider the following desingularized Laplace transform 

c/)(t,s)= [ {e~ sx - 1 + sx)u t (dx). (2.20) 
Je 

We substitute / = e~ sx — 1 + sx in the equation of evolution (2.5) to find 

d t ^={f\L{u t ))= l -{d s ^) 2 . (2.21) 

Equation (2.21) is the Hamilton-Jacobi equation associated to the inviscid 
Burgers equation. Thus, we let 

iP{t,s)=d s <P= [ (l-e- sx )xv t (dx). (2.22) 
Je 

Then from (2.21) we have 

d t i/j - ipd s ip = 0. (2.23) 

The exact solution to (2.23) with initial data tpo(s) may be found by the 
method of characteristics. The characteristic originating at sq is denoted 



s(t, s ) = s - ipo(s )t. 
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Let t(so,si) denote the time for two characteristics originating at so < si 
to intersect. Then, if this is the first intersection 

1 Vo(si) - ^o(so) , 1 1 

whence < r(so,si) < 



i s\-s Q d s ip (s) ' ^V'o(si)' 

where the inequalities follows from the mean value theorem and the complete 
monotonicity of d s ipQ. Thus, letting sq = and s\ — > 0, we see that the least 
time taken for characteristics to intersect is given by 

T g J = d s M0) = [ x 2 u (dx) = m 2 (0). 

J E 

Without loss of generality, we may assume that the initial data is normalized 
so that m2(0) = 1 = Tg e i and thus < t < 1. This normalization assumption 
is analogous to the assumption that mi = 1 for the additive kernel. Equation 
(2.23) should be compared with equation (2.3). In fact, given initial data 
V'o('S), by changing the time scale in (2.3) it is easy to check that ip(t,s) is 
the solution to (2.23) if and only if 

^(t,s) = Y3^(-log(l -*)>*)> (2.24) 

where ip(t,s) is the unique solution to (2.3) with initial data ipo. The next 
lemma follows immediately from Lemma 2.9. 

Lemma 2.10. Suppose vq(x) G M.+ with f E x 2 vo(dx) = 1 Then equation 
(2.23) determines a map [0, 1) 3 1 1— > i/ t g A4 + smc/i £/iat 

L m 2 (t) = j E x 2 u t {dx) = (1-t)" 1 . 

x 2 fj is weakly continuous on [0,1). 

It is natural to term the measure zy t the solution to Smoluchowski's co- 
agulation equation with kernel K = xy. However, it is harder to formulate 
a completely natural well-posedness theory in this case, and we will settle 
for a reasonable compromise. 

Definition 2.11. Define the norm 

sup + := n/ii^ (2 . 2 5) 

x,y>o xy 



and the associated Banach space V = {/ G Co(£?)|||/||v < °°}- 
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It is clear that V is a Banach space. The norm || • ||y is natural in the 
following sense. 

Lemma 2.12. Let v G M.+ such that rri2 < oo, and let L(v) be defined by 
(2.5) with K = xy. Then L(v) defines a bounded linear functional on V 
with norm < m\/2. 

Proof. Since \ f(x + y) - f{x) - f(y)\ < \\f\\ v xy we have 

\(f,L(u))\<U f \\f\\ v xW^dx)v(dy) = 7 ^\\f\\ v . 

J E J E 

□ 

It is easy to check that finite sums f(x) = J2k=i a h(l ~ e~ SkX ) are in V. 
We would like to believe that these functions are dense in V, but this seems 
hard to prove, as the norm above is unwieldy. Instead we will work with C 2 
functions and use the following, whose easy proof we omit. 

Lemma 2.13. Let f be aC 1 function such that f(0) = and f is Lipschitz. 
Then\\f\\ v <2Up(f). 

Theorem 2.14. Let v G M + satisfy m 2 (0) < oo. Then Smoluchowski's 
coagulation equation with kernel K = xy has a unique solution with initial 
data vq on the time interval [0, m2(0)~ 1 ). 

Proof. Without loss of generality we may suppose that m2(0) = 1. The 
measures v t of Lemma 2.10 are a candidate solution, and it is easy to check 
that conditions (1) and (2) of Definition 2.3 are satisfied. Since m2(t) = 
(1 — t)' 1 by Lemma 2.12 one sees that L{v t ) is a bounded linear operator 
on V and in particular 

(l-8x-e- ax ,L{ Vt )) = ~{</> 8 ) 2 . 

Thus, <p(t,s) solves (2.21). Yet, some care is needed in checking that (2.8) 
holds in full generality. Let / £ f 2 . We apply Lemma 2.2 to /' (notice, 
not f) to obtain an approximation P n (x;s) with sup^, \P n — f'\ < e, and 
sup^ l-P/j — f"\ < e. Observe that, we may rewrite 

n n 

P n (x; s) = J2 On,*(l - e- sx ) k = E ^M 1 ~ e ~ SkX ) 

k=l k=l 
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by expanding (1 — e ~ sx ) k with the binomial formula, and defining b n ^ as 
the corresponding linear combinations of a n ^. We integrate P n to obtain 
the approximation 

n , 

k=i S 

Notice that P n (0; s) = 0. Therefore, by the fundamental theorem of calculus, 
we also have 

\P n (x; s) - f{x)\ < f \P' n {z- s) - f"(z)\dz < ex, 
Jo 

and upon integration again, 

sx 2 

\f(x)-Q n (x;s)\ < — . 

But we then have, 

\{f-Q n ,v T )\ < ^J^x 2 v T (dx) = £ -m 2 {r), r G [0,1). 

Since sup^, \ f" — Q^\ < e we apply Lemma 2.12 and Lemma 2.13 to obtain 

\(f-Q n ,L(u T ))\<em 2 2 (T), re [0,1). 

This proves that v t is a solution. 

It is slightly harder to prove uniqueness in this case. As in Theorem 2.4 
and Theorem 2.8 it suffices to deduce uniqueness of the measure valued 
solution via uniqueness of solutions to (2.21). The obtstruction is that 
it is not clear from the definition of the weak solution that (2.21) holds, 
since the test functions 1 — sx — e~ sx do not lie in £ 2 . This can be over- 
come with an approximation argument that we only sketch. We consider 
f n (x) = (1 — sx — e~ sx )xn(x) where \n is a C°° cut-off function such that 
Xn = 1) x < n, Xn = 0) x > n + 1. By the monotone convergence theorem, 
lim n _ >00 (/ n , Vt) = (1 — sx — e~ sx ,ut). By Lemma 2.12 and Lemma 2.13 
{f n ,L(v t )) is well-defined and uniformly bounded by C s rri2(t) 2 . We may 
then use the dominated convergence theorem to deduce that (2.21) holds in 
the limit n — > oo. Uniqueness of v t follows. □ 
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3 Regular variation 

Several formal calculations by physicists working on Smoluchowski's equa- 
tions take the following form : (1) assume that the number density n(x) ~ x a 
for some scaling exponent a, and (2) conclude based on physical arguments 
that a takes a particular value. The theory of regular variation helps us 
makes these formal calculations precise, and lays bare the mechanism con- 
trolling the approach to scaling form. Our primary source is Feller's book, 
and we restate below useful results from [16, VIII. 8]. The theory of reg- 
ular variation has many applications in analysis and probability, and an 
authoritative text, rich in examples, is [9]. 



3.1 Rigidity of scaling limits 

Loosely speaking, a function is slowly varying if it is asymptotically flat 
under changes of scale. Precisely, we say that a positive function L(x) is 
slowly varying at infinity if 

lim 4r^r = 1> for a11 f > °- (3- 1 ) 

x-^oo L[x) 

For example, all powers and iterates of log x are slowly varying at infinity. 
If we consider the limit x — > instead, we obtain functions that are slowly 
varying at 0. 

A function N{x) is regularly varying at infinity with index p € K if there 
is a slowly varying function L(x) such that 

N(x) ~ x p L{x) as x -> oo. (3.2) 

The notation ~ means lim^oo N(x)/x p L(x) = 1. 

The notion of regular variation is intimately related to necessary and 
sufficient condtions for the existence of scaling limits. This is reflected in 
the following classical "rigidity" lemma [16, Lemma VIII. 8. 3], which will be 
one of our principal tools. 

Lemma 3.1. Suppose that a n+ \/a n — > 1 and X n — > oo as n — > oo. If ip is a 
positive, monotone function such that 

lim a n ip ( — I = g(s) < oo 

exists for s in a dense subset of (0, oo), and g is finite and positive on some 
interval, then ip varies regularly at and g(s) = cs p for — oo < p < oo 
some c > 0. 
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3.2 Tauberian theorems 

We will rigorously deduce the asymptotics of v by the beautiful Hardy- 
Littlewood-Karamata Tauberian theorem [16, XIII. 5]. 

Theorem 3.2. If L is slowly varying at infinity and < a < oo, then the 
following are equivalent: 

v(x) ~ x a L(x) as x — > oo, 

and 

n(s) ~ s~ a L Q T(l + a) as s -► 0. 

Moreover, this equivalence remains true when we interchange the roles of 
the origin and infinity, namely when s — > oo and x — > 0. 

We will use the following lemma to show that there is no loss of generality 
in working with ip instead of n. 

Lemma 3.3. Suppose d s ip is the Laplace transform of a positive measure. 
Let a < 1 and L be a function slowly varying at 0. The following are 
equivalent. 

1. - ^(0) ~ s l ~ a L(s) as s ^ 0. 

2. d s ip{s) ~ (1 - a)s- a L(s) as s -► 0. 

Proof. Suppose (1). Since ip(s) — ip(0) = s a L(s)h(s) with lim s _>o h(s) = 1, 
without loss of generality we may write ip(s) — ip(0) = s a L(s). Fix a > 1. 
Then by the mean value theorem and the complete monotonicity of d s ip we 
have 

s(a - l)d s 4>(s) > V(as) - ^(s) = s^L^s) (a 1 ^ 1 ^^- - 1 



L(s) 



Thus, letting s — > 0, and using (3.1) we have 



s^o L(s) a — 1 

Since a > 1 is arbitrary, we may maximize the right hand side to obtain 

liminf - > 1 — a. 

s->o L(s) 
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Choosing a < 1 and using a similar argument yields, 

hm sup — — < 1 — a. 

s ^o L(s) 

Thus, a s ^(s) ~ (1 - a)s- a L(s). 

Conversely, assume (2). Then we have 

^(s) - ^(0) = (1 - a) /" t- a L(t)dt = (1 - q)s 1 - q L(s) dt - 

JO JO M s ) 

Since L is slowly varying at zero, then for any constants A > 1 and <5 > 
there exists so such that for < s < sq, < t < 1 we have L(st)/L(s) < 
At~ 6 (this is not hard to show, but see Theorem 1.5.6 in [9]). Then (1) 
follows by dominated convergence. □ 



4 Scaling solutions for the constant kernel 

4.1 Mittag-Leffler distributions 

The scaling solution 

n(t, x) = t~ 2 exp(-z/t), t > (4.1) 

is the continuous limit of a special solution found by Smoluchowski [4]. Kreer 
and Penrose proved that the rescaled number density t 2 n(t, xt) converges 
uniformly to e~ ax on compact sets, under the assumption that the initial 
number density no (a;) be C 2 and have exponential decay in x [22]. The 
constant a is determined by the initial mass. 

In this section we show that the solution (4.1) is just one of a one- 
parameter family of scaling solutions given by 

n(t,x) = t- 1 - 1 / p n p (xr 1 / p ), t > 0, p e (0, 1], (4.2) 

where n p (x) = F' p {x) is the density, and F p the distribution function for the 
Mittag-Leffler distribution 

Of these solutions, only the solution (4.1) for p = 1 has finite mass, and 
the others have fat tails. The Mittag-Leffler distribution was studied by 
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Pillai [31] who showed that these distributions are infinitely divisible and 
geometrically infinitely divisible for p £ (0,1]. 

For our purposes, it is especially relevant that the Mittag-Leffler distri- 
bution has Laplace transform 



/'OO /'OO -1 

/ e- sx n p (x)dx= / e - sx F p {dx) = -—. 
Jo Jo 1 + s 

In terms of the Mittag-Leffler function 



(4.4) 



one has F p (x) = 1 — E p {—x p ). It is interesting and useful to note the 
following recent calculation of Tsoukatos [33] . 

Lemma 4.1. For < p < 1 we have 

71 Jo v + cos it p) z + (sm7rpj z 

Hence E p (—x p ) and n p (x) = —d x E p (—x p ) are completely monotone. 
Proof. We sketch the argument of Tsoukatos [33]. Since (4.4) implies 



f 

Jo 



one can invoke the Laplace inversion formula and evaluate it by deforming 
the contour to fold along the negative real axis to obtain, for any a > 0, 



,a+ioc s p ds U f°° - rr r p e inp dr 

e 1 



E p (-x p ) = / e M -i— - = -Im / 



1 _)_ r p e mp r 

and the result follows. □ 

Note that the complete monotonicity of E p (—x) (conjectured by Feller 
and proved by Pollard in 1948) also implies the complete monotonicity 
above, since x p is positive with completely monotone derivative. A point of 
confusion in the literature is that the term "Mittag-Leffler law" is used by 
some for the distribution whose Laplace transform is E p (— s) [9]. 
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4.2 Scaling solutions 

Let us now check that n p (t, x) defined by (4.2) are indeed solutions to (1.1) 
when K = 2. We take the Laplace transform of (1.1), and its limit at s = 0, 
to obtain 

§2 =,•-!*«, OH ^ = -,M)2. (4.7) 

In analogy with (4.1) we make the ansatz rj(t,s) = t~ 1 ri p (s\(t)), r/ p (0) = 1 
in equation (4.7). Letting £ = sA we have 

(t^tf P = -ri P 0--Vp)- (4-8) 

Equation(4.8) may be simplified by separating variables. We let p denote 
the separation constant to obtain 



^ = -pn p {\ - T, p ), j = y t (4.9) 



The general solution to (4.9) is 



^) = TT-5;' A(t) = (c 2 t)VP (4.10) 

where Ci, c 2 > are arbitrary constants. We combine the two solutions to 
find for each p, a family of solutions related by scaling in s and t, 

" (t ' s) = t-1 TT^«' f>0 ' 

rj(t, s) is completely monotone if and only if p £ (0, 1] [31], thus it is only for 
p £ (0, 1] that we obtain positive solutions to equation (1.1). By a trivial 
scaling we may achieve c\ = C2 = 1, and then n(t,x) is given by (4.2). The 
scaling solutions have finite mass only when p = 1, and in this case the mass 
is conserved. 

4.3 Asymptotics of scaling solutions 

We may use equation (4.3) to obtain the convergent expansion 

~ (_l)k+l pk-l 
W ' (g) = g - T(pk) ' X> ° (412) 
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which implies the divergence n p (x) ~ x p ~ 1 /F(p) as x — > 0+, for p G (0, 1). 

Since n p (x) is completely monotone for p £ (0, 1), its asymptotic prop- 
erties as x — > oo may be obtained rigorously by differentiating the formula 
in Lemma 4.1 and using the Tauberian theorem. We obtain 

n p {x)~x-v- l T{2 + p)^^ = ^-^- asx^oo, (4.13) 

?r(i + p) -r(-p) 

using y = 1 + p in the identity 

r(l+y)r(l-y)^^ = l. (4.14) 
ny 

From (4.13), or because — n p (s) ~ ps p_1 as s — > 0, we find 
/" x px 1_/9 

y yn p (y)dy ~ r ^ 2 _ ^ » as x ^ oo. (4.15) 

The case p = 1/2 curiously admits the exact solution (see [1, Ch. 29]) 

n /a( x ) = ; ~~ e^erfc-v/x. 

\/7rx 



5 Weak convergence for the constant kernel 

Let vi be the measure-valued solution of Smoluchowski's equation with ker- 
nel K = 2 obtained in section 2.2, given initial size distribution vq that is 
a finite measure. We normalize to define a probability distribution function 
(the size-biased distribution) 

n*,*) = i#% *>0. (5.1) 
V ; ^((0,oo))' V ; 

We then have the following characterization of permissible limits under 
rescaling and their domains of attraction. Below, we call a probability dis- 
tribution function F*(x) nontrivial if F*(x) < 1 for some x > 0, meaning 
the distribution is proper (lmx^oo F(x) = 1) and not concentrated at zero. 

Theorem 5.1. 1. Suppose there is a rescaling function \{t) — ► oo and a 
nontrivial probability distribution function F*{x) such that 



lim F(t,X(t)x) = F*(x) 

t—>oo 



(5.2) 
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at all points of continuity of F*. Then, there exists p G (0,1] and a 
function L slowly varying at infinity such that 

rx 

I yvo{dy) ~ x 1 ~ p L(x) as x — > oo. (5-3) 
Jo 

2. Conversely, suppose there exists p £ (0, 1] and a function L slowly 
varying at infinity such that (5.3) holds. Then it follows that there is 
a strictly increasing rescaling A(t) — > oo such that 

lim F(t, X(t)x) = F p {x), x G (0, oo), 



t— >oo 



where F p is the Mittag-Leffler distribution function defined in (4-3). 

Proof. Our proof is based on the Laplace transform. We first reformulate 
(5.2) and (5.3) in terms of (po(s). Firstly, by the well-known characterization 
of weak convergence by the Laplace transform [16], (5.2) is equivalent to the 
assertion that the Laplace transforms converge pointwise, i.e., that 

lim ^' sA ~ 1) ^ r,*(s) := f°° e'^F^dx), s G [0, oo). (5.4) 

The assumption F*(x) < 1 for some x > ensures that < 77* (s) < 1 for all 
s > 0. Since n(t, s) = <p(t, 00) — (p(t, s) and <p(t, 00) = n(t, 0), by the solution 
formula (2.9) we have 

vit^X' 1 ) = 1 + yo(5A~ 1 )y?o(oo)~ 1 
i?(t,0) l + t^o^A- 1 ) ' { -> 

Because <po(0) = 0, existence of the limit in (5.4) with < 77* (s) < 1 is 
equivalent to the existence of 

g(s) := lim ^(aAC*)" 1 ) (5.6) 

t^oo 

with < g(s) < 00, for all s > 0. 

The behavior in (5.3) may also be reformulated in terms of the Laplace 
transform. Applying Theorem 3.2 and Lemma 3.3 to the mass distribution 
function appearing on the left-hand side, we find (5.3) equivalent to 

Ms) ~ r(2 ~ P) s p £(-), as s 0. (5.7) 
p s 

Now to prove the first part of the theorem, we take t = 1, 2, ... in (5.6) 
and apply Lemma 3.1 to conclude that <po(s) is regularly varying at and 
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g(s) = cs p for some c > and p > 0. In fact /? > since r/*(s) = (1 + cs p ) _1 
must satisfy r/*(0) = 1. Hence (5.7) holds, and it remains to show that 
p £ (0, 1]. This will follow from complete monotonicity of the limit. Since 
rj(t,s) is completely monotone, it follows that r/*(s) is completely monotone 
since it is the limit of a sequence of completely monotone functions. This is 
possible only if p G (0, 1], since the second derivative of (1 + cs p )~ l is not 
ultimately positive if p > 1 [31]. 

Conversely, we prove the second part by showing that (5.7) implies (5.4) 
with = F p . We define X(t) for sufficiently large t by 

tMKty 1 ) = i- (5-8) 

X(t) is strictly increasing because (fo(s) is strictly increasing. Moreover, 
limt^oo X(t) = oo since ip(0) = 0. Since (fo is regularly varying with index p 
we have 

lim t m (s\- r ) = lim ^ o( / S x A "' ) = s p . 

t^oo t-^oc (po[X ) 

But then (5.5) yields (5.4) with 77* (s) = (1 + s p ) _1 . □ 

Remark 5.2. Let ipo(s) = s p L{l/s). Then equation (5.8) shows that 

X(t)L(X(t)) 1/p = t 1/p . 

Comparison with the time scaling X{t) = t 1 ^ for the self-similar solution 
(4.10) shows that X(t) chosen in the proof is essentially the time scaling of 
the self-similar solution, possibly modified by a slowly varying correction. 

Remark 5.3. When p = 1 the condition for being attracted to the expo- 
nential distribution is f£ yvo{dy) ~ L(x) as x — > 00. Thus, all solutions with 
initially finite mass are attracted to the finite-mass exponential distribution, 
but it is not necessary for this that the initial mass be finite. It suffices that 
the mass distribution function diverge sufficiently weakly. 

Remark 5.4. A remaining nontrivial possibility to discuss is that a nonzero 
in (5.2) may exist where the function is a defective probability distribu- 
tion, satisfying F*(oo) < 1. If this is the case, then since ??*(0 + ) = F*(oo) < 
1 it follows g(s) = cs p with p = 0, and that <po(s) ~ L(l/s) is slowly 
varying at 0. We cannot ensure (5.3) in this case. Instead we note that 
ifo(s)/s = / °° e~ sx vo(dy) dx, and it follows from the Tauberian theo- 
rem and the fact that x i-> v${dy) is monotone ([16], XIII. 5. 4) that the 
tail distribution function is slowly varying at 00, with 

/>oo 

v ((x,oo)) = / vo(dx) ~ L{x). (5.9) 

J x 



28 



Conversely, if (5.9) holds, then tpo(s) is a function slowly varying at that 
strictly increases. For any c £ (0, oo), we can choose A(i) strictly increasing 
such that t(^(A(t)~ 1 ) = c. Then it follows that (5.4) holds with r/*(s) = 
(1 + c) _1 for s > 0, so (5.2) holds with the defective distribution function 
F*(x) = (l+c) _1 . This means that under such scalings, an arbitrary fraction 
of the particle sizes concentrate at and the rest escapes to infinity. 

6 Scaling solutions for the additive kernel 
6.1 A one-parameter family of solutions 

Golovin found an exact solution to Smoluchowski's equations with monodis- 
perse initial condition for K = x + y. One may take limits in his solution to 
obtain the scaling solution [4] 



This solution has sometimes been criticized as unphysical, since the number 
of clusters is infinite. However, recently Deaconu and Tanre proved a result 
equivalent to weak convergence to this solution, under restrictive assump- 
tions on the initial data (the existence of all moments and their domination 
by the moments of a Gaussian random variable) [12]. 

We will consider only solutions of finite mass, normalized to 1. The 
solution (6.1) is but one of a one-parameter family of solutions (independent 
of the trivial scaling c 2 n(t, cx)). For each p 6 (0, 1], in this section we derive 
finite-mass scaling solutions in the following form, with (3 = p/(l + p): 



n(t, x) 



1 



,x" 3/ V exp(-e" 2 *x/2) 



(6.1) 



V2^ 



n(t,x) =e- 2t/f3 n p (e- t ^x), 



(6.2) 




(6.3) 



k=l 



The associated mass distribution function is given by 



M, 



r P (x) = / yn p (y) dy = Y^ 



(-l)^- 1 ^ 
k\ 



r(l + k — kf3) 



sin nk(3 



nk(3 



(6.4) 



Remark 6.1. It is an interesting fact that these scaling solutions are related 
by a nonlinear scaling to the Levy stable laws in probability theory. Feller 
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([16], XVII.7) gives the formula 

( JvTT 

p(x; a, 7) = — -JT T ^ + k,a) Sin 2^ (7 " a) (6 ' 5) 
k=i 

for a family of stable densities, for 1 < a < 2, \j\ < 2 — a. Taking a = 
(1 — 7 = 2 — a, we find that the mass density from (6.3) satisfies 

xn p (x) = x l3 ~ 1 p(x lS ; 1 + p, 1 — p). (6.6) 

These remarkable self-similar solutions were first discovered by Bertoin [7] 
(then independently by us). Bertoin's derivation explains the nonlinear 
rescaling formula in terms of a scaling property of Levy stable processes, 
and he writes the self-similar solution in the form 

n(t, x) = e'V-y e'V; 1 + p, 1 - p). (6.7) 

It is quite remarkable that there are two scaling limits associated to this 
solution. One of them is more transparent in (6.7): we find that e t n(t,x)dx 
converges vaguely towards the measure x^~ 2 dx. This is the scaling limit 
alluded to in Corollary 1 of Bertoin's article [7], and it suffices to uniquely 
identify the self-similar solution in the class of eternal solutions. On the 
other hand, (6.2), reflects more clearly the self-similar nature of the solution 
relative to the mean cluster size e 1 ^ . 

Note that the stable densities are defined on the whole line (—00,00). 
We obtain total mass 1 on (0, 00) through the nonlinear rescaling. Moreover, 
if F(x; a, 7) denotes the distribution function for the stable law with density 
p(x; 0,7), then the tail of the mass distribution corresponds to this through 

1 - M p (x) = 0~\l - F( X P; l + p,l- p)). (6.8) 

The total number of clusters diverges for all the solutions in (6.2). This 
is caused by the predominance of small clusters, and it may be desingular- 
ized by working with the variable ip introduced in Section 2.1. The scaling 
solutions above are given implicitly in terms of <p p = tp p {s) satisfying 

* = + (6.9) 

Since u p (s) = d s ip p is the Laplace transform of the mass distribution, differ- 
entiating (6.9) we have 

1 
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which exhibits a connection to the Mittag-Leffler distribution (see (4.4)). 
We use this to show below that the mass distribution is infinitely divisible. 
When p = 1, equation (6.9) is a quadratic equation with two solutions, one 
of which is <p(s) = y/l + As — 1, which corresponds to (6.1). For p G (0, 1) we 
will solve (6.9) by Laplace's inversion formula as an infinite series, to obtain 
(6.3). 

6.2 Scaling solutions 

We may derive mass-preserving scaling solutions to (2.3) as follows. Let X(t) 
be a rescaling to be determined, and let £ = sX. We substitute the ansatz 
f(t,s) = X^ 1 ip p (sX) = A~V P (£), in (2.3) to obtain 

-j(<Pp - Zfy<P P ) - <Ppfy<P P = -<P P - (6.11) 

We separate variables by letting A/A = a or A = c\e at . Then by (6.11) 

« - <Pp)dz<p p + (1 - a)ip p = 0. (6.12) 

Equation (6.12) is not separable, but it may be solved implicitly by rewriting 
it as the linear equation 

* a C 1 (6.13) 



difp a — 1 ip p 1 — a 

Put p = (a — l)" 1 so a = (1 + p)l p = 1/(3. Integrating, we find a family of 
nontrivial solutions determined by 

£ = ¥> P + C2¥>£ +p , c 2 >0. (6.14) 

The range of admissible p is narrowed by requiring that lim^o = 1 
(finite mass) which implies p > 0. Without loss of generality we may take 
C2 = 1, since we can recover all other solutions by a trivial scaling. 

We now show that p > 1 is inadmissible. Let £/ (£) = </?(£)/£. Integrating 
(2.1) by parts we see that U(£) = J °° e~^ x N(x)dx where N(x) = v((x, oo)) is 
the tail distribution. In particular, £7(0) = 1 and [/ is completely monotone. 
Dividing (6.14) by £ and differentiating we see that 

U'(£) = ^ > as £ -f 0, 

which is impossible if J7 is completely monotone. Thus, the admissible range 
of non-trivial solutions is restricted to p 6 (0, 1]. 
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6.3 Series expansion and asymptotics of scaling solutions 

The asymptotic properties of the scaling solutions for p G (0, 1) may be 
rigorously obtained from Theorem 3.2. By (6.9), (p" p ~ p(p + l)s p_1 as 
s — ► 0. But, ip'p = J °° e~ sx x 2 n p (x)dx, and Theorem 3.2 implies that 

J y 2 n p {y)dy ~ Y(2 + - l p) xl ~ P ' as x ^ °°- ( 6 - 15 ) 

Thus the second moment is finite only for p = 1. For < p < 1 the mass 
distribution has fat tails. Equation (6.15) is a weak version of the pointwise 
behavior 

n p {x) ~ |r( + ^| x ~ (2+p) asa^oo, p G (0,1), (6.16) 

which follows from (6.6) due to the known power-law asymptotics of the 
stable densities [5]. 

The behavior as x — ► is described completely by the series (6.3), derived 
as follows. We rewrite (6.9) in terms of U = <p p /s and (3 = p/(l + p) as 

U{s) = s-P(l - U) 1 -! 3 . 

We solve for U using Lagrange's inversion formula (see e.g. [20, 6.3] for a 
similar calculation), obtaining 

00 -k@ fik-l 

U ^ = E^Tdx^( F ^ k \^ with F(x) = (l-x)^. 

k=l 

We evaluate the derivatives and find that 

oo — fc/3 k 

d^ p (s) = d s (su) = ^2 V (_i)fe_1 no - ( e - 17 ) 

fc=l ' j=l 

This is the Laplace transform of the mass distribution function given through 
term-by-term Laplace inversion as 

k=l y 1 ' j=l 

We then deduce (6.4) using y = k[3 in the gamma-function identity (4.14). 
By differentiating (6.4) we obtain the number density in (6.3). 

It is straightforward to check that when p = 1, the even terms vanish and 
(6.3) reduces to the function (47r)~ 1 / 2 x~ 3 / 2 e _:r//4 which is a scaled version 
of (6.1). Correspondingly, M\(x) = erf(iy / x). One may also check that the 
series solution above is absolutely convergent for x £ (0, oo). Thus, n p is 
analytic. 
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7 Weak convergence for the additive kernel 

We let v t be a solution of Smoluchowski's equation with kernel K = x + y 
given initial data vq with finite mass normalized to / °° x fo(dx) = f . Then at 
all times the mass distribution is a probability distribution, with distribution 
function denoted 

M(t,x)= [ yv t {dy). 
Jo 

It follows from the explicit solution (2.16) that the Laplace transform of the 
mass distribution satisfies u(t, 0) = 1 for all t > 0, and lim^oo u(t, s) = for 
s > 0. This phenomenon of concentration is equivalent to the assertion that 
asymptotically all the mass escapes to infinity. As earlier, we may hope that 
suitable rescaling in s will give convergence to a nontrivial limit. Precisely, 
we have the following characterization. 

Theorem 7.1. 1. Suppose there is a rescaling function \(t) — ► oo as 
t — > oo and a nontrivial probability distribution function M*(x) such 
that 

lim M (t, \{t)x) = M*(x) (7.1) 

at all points of continuity of M*. Then there exists p G (0,1] anc? a 
function L slowly varying at infinity such that 

/ y 2 uo(dy) ~ x 1_p L(x) as x — > oo. (7.2) 
J o 

Conversely, assume that there exists p £ (0, 1] and a function L slowly 
varying at infinity such that (7.2) holds. Then there is a strictly in- 
creasing rescaling X(t) — > oo smc/i i/ia£ 

lim M(t, X(t)x) = MJx), < x < oo, 

where M p from (6.4) is the mass distribution function for a scaling 
solution. 

Proof. We will prove the theorem after reformulating (1) and (2) as equiv- 
alent assertions using the Laplace transform. Firstly, the weak convergence 
of the mass distribution M(t, \x) is equivalent to the pointwise convergence 
of its Laplace transform 

poo 

lim u(t,s\- 1 ) = u*(s) := / e~ sx M*(dx), < s < 00. (7.3) 
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The assumption M*(x) < 1 for some x > ensures that < u*(s) < 1 for 
s > 0. Secondly, (7.2) is equivalent to — d s u ~ s f) ~ 1 L(l/s)T(2 — p) as s — > 
by Theorem 3.2. Since p € (0, 1], by Lemma 3.3 this is equivalent to 

1 - u {s) ~ s^L f-^j F(2 - p) as s 0. (7.4) 



We prove the first part of the theorem by showing that (7.3) implies (7.4). 
Since u*(s) is a limit of completely monotone functions, it is itself completely 
monotone. Moreover, since u(t, s) = d s (p we also have the convergence 

lim \ip(t, sA" 1 ) = </?*(s) = / u*(s')ds'. (7.5) 
t->oo ,/ 

Clearly, 99* is strictly increasing. 

In what follows, we consider a(t, s) defined by replacing s with s/X(t) in 
(2.15), i.e., by 

a-ip (a) = sX- 1 - (p(t,s\- L ). (7.6) 
It then follows from (7.6) that as t — > 00 with s fixed, we have a — > 0, and 

lim A(cr — </? (<t)) = s — </5*(s). (7.7) 

From (2.14) and (7.5) we then have 

lim Ae _ *^>o(c) = lim Xip(t, sA _1 ) = <£>*(s). (7-8) 

Replacing s by sA _1 in the exact solution (2.16) we also have 

lim — — ^ — = uJs). 

t^oo e *(l - u (<r)) +u (a) 

Since a — > and uo(0) = 1 we deduce that 

lim e t (l-u (a))= 1 ~ U *[ s) . (7.9) 

We now show that a — > at the rate a(i) := e*A(t) _1 (lim^oo a(i) = by 
(7.8)). We may rewrite (7.6) as 

a = a(t) [(p*(s) + se~ l + Ae _ Vo(o")(l - e~*) - ip*(s)] 
= a(t)[ip*{s) + r(t,s)], 
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where the error term r(t,s) — > by (7.8). Therefore, a is asymptotically a 
scaling of </?*(s). We now claim that for all s > 0, 

lim e*(l - uoKt)^))) = (7.10) 

t^oo U*(S) 

Fix s > 0, and let 5 > be sufficiently small. Since <p* is strictly increasing, 
for sufficiently large t (depending on s and 5) we have 

</?*(s — S) + r(i, s - 5) < <p*(s) < tp*(s + (5) + r(i, s + 5), 

whence 

1 — uo(a(t, s — 5)) < 1 — uo(a(t)<p*) < 1 — uo(a(t, s + (5)). 

Multiply by e* and take i — ► oo, then <5 — > 0. The claim (7.10) then follows 
from (7.9). 

It follows directly from (7.10) and Lemma 3.1 that 1 — uq is regularly 
varying at with some exponent p G R, and the limit in (7.10) has the form 
cipZ for some positive constant c. Clearly < p < 1 since uo is bounded 
and completely monotone and u*(0) = 1 by the hypothesis that M* is a 
probability distribution. This finishes the proof of the first part. Note 
furthermore that (6.10) holds after scaling s. 

We prove the converse statement by showing that (7.4) implies (7.3) with 
u* = Up. By the explicit solution formula (2.16), it suffices to show that as 
t — > oo we have 

e'(l - u (a(t, sX' 1 ))) = u(t, sX' 1 )- 1 - 1 (1 + p)^(s) p (7.11) 

for all s > 0. We write a(i) = e*A(t) _1 as earlier. We choose X(t) to satisfy 
A(0) = and 

e t (l-« (a(t))) = l + p, t>0. (7.12) 

Then a(t) — > 0, A(i) is strictly increasing and A(i) — > oo. It follows from 
(7.4) that for fixed ip* > 0, as t — > oo we have 

e*(l - ««,(»(%,)) = (1 + />) V^Tr + ^ ^ 

1 — uo{a{t)) 

For fixed 99* > we define s(t,ip m ) as the value of s determined from (7.6) 
using a = a(t)ip*. Then it follows that for all ip* > 0, 

lim u(t, s{t, ip^X- 1 )- 1 - 1 = (1 + p)ipP. (7.14) 



35 



Using (2.14) with (7.6) we have 

s(t,(p*) = X{aip* — <po(a<p*)) + Ae _ Vo(a<^*) 

= (1 + p)— t: rrr- + "WW*. 

a(l — -uo(aj) a 

Since 1 — uo ~ s p L(l/s), the proof of Lemma 3.3 shows that s — po(s) ~ 
(1 + p)~ 1 s p+1 L(l/s). Therefore, as t —>■ oo we have 

»(«. ~ ^czy + y - + (7as) 

Now to prove (7.11), fix so > 0. Then there is a unique ip* = v? p (so) > so 
that so = i P*( i P* + !)• Since ip* 1— > s(t, (/?*) is strictly increasing in 99* for all 
i, by substituting tp* ± 5 for in (7.14) we easily deduce (7.11). □ 

Remark 7.2. When 1 — uo = s p , the choice of time scaling X(t) in (7.12) 
gives X(t) = e( 1+1 / p )* in accordance with (6.12). More generally, when 
1 — no = s p L(l/s) the rescaling X(t) is modified by a slowly varying correc- 
tion. The choice of time scale when p = 1 and the second moment is finite 
deserves special comment. In this case we find X(t) = e 2t . In the applied lit- 
erature it is common to define mean cluster size as a ratio of moments. Let 
mk(t) = Jq°° x k i>t{dx). It is clear that any ratio of the form m^+iit) / 'mk(t) 
has the dimensions of length, and two distinct, but natural definitions of 
mean cluster size are (see [34]) 

cut) = rr and C2U) = --. 

U 7 mo(t) K ' mi(i) 

For Golovin's solution the initial data are monodisperse, and mo(0) = 
mi(0) = m2(0) = 1. It is easy to calculate explicitly that c\{t) = e t , 
but C2(t) = e 2t . Thus the two notions of cluster size differ, and only e 2t is 
the correct scaling for convergence to self-similar form. More generally, e 2t 
is the only choice of time scaling that fixes both mi and 7772 as required for 
convergence to self-similar form. 

Remark 7.3. Theorem 7.1 with p = 1 shows that the mass distribution 
is attracted to the classical self-similar solution for all initial data vq with 
finite first and second moment. But for this behavior it is not necessary 
for the initial data to have finite second moment. It suffices that it diverge 
sufficiently weakly, with J* y 2 v$(dy) ~ L(x) slowly varying at infinity. 



36 



Remark 7.4. A remaining nontrivial possibility is that a nonzero limit in 
(7.1) may exist with M* defective, satisfying M*(oo) < 1. If this is true, 
then most of the proof of the first part of the theorem carries through. The 
limit in (7.10) must have the form c<p* with c > 0, but we must have p = 0, 
since u*(0) < 1. Moreover it follows 1 — uq(s) ~ L(l/s) is slowly varying. 
We do not obtain (7.2) in this case. Instead, we note 



uq(s) 



f'OO /'OO 

/ e~ ax / vMdy)dx. (7.16) 

JO Jx 



As for the constant kernel, it follows from the Tauberian theorem and mono- 
tonicity that the tail of the mass distribution is slowly varying at infinity, 
with 

yVQ {dy)~L{x). (7.17) 



f 

J X 



In the converse direction, if (7.17) holds, then 1—uq(s) is a strictly increasing 
function that is slowly varying at 0. For any c £ (0, oo) we can choose X(t) 
strictly increasing so that e*(l — uo(a(t))) = c. Then it follows as in the proof 
of the second part of the theorem that (7.3) holds with u*(s) = (1 + c) _1 
for s > 0, so (7.1) holds with the defective distribution function M*{x) = 
(1 + cT 1 . This means that under such scalings, an arbitrary fraction of the 
mass concentrates at and the rest escapes to infinity. 

We conclude this section with a useful observation about the self-similar 
solutions. 

Theorem 7.5. For each p G (0,1], the probability distribution M p is in- 
finitely divisible. 

Proof. It suffices to show that the Laplace transform u p = e~^ p where 
VV>(0) = 0, and ijj p has completely monotone derivative [16, XIII. 7.1]. By 
(6.10), V P = log(l + (1 + p)<p$). Clearly, ip p (0) = 0. Moreover, 

(1 + p) P ip p f r 1 u p 



i + (i + pH 



By Theorem 7.1 u p is completely monotone. The other factor can be written 
as a composition with the Mittag-Leffier distribution 

The function s p ^ 1 is completely monotone, as is (1 + s p ) _1 . Thus, their 
product is completely monotone. Thus, the composed function above is 
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completely monotone since it is the composition of a completely monotone 
function with a function that has a completely monotone derivative [16, 
XIII. 4. 2]. Finally, d s ifj p is the product of two completely monotone functions, 
and is hence completely monotone. □ 



8 Approach to self-similar gelation for the multi- 
plicative kernel 

8.1 McLeod's solution 

McLeod found the following explicit solution to the discrete Smoluchowski 
equation (1.1) for K = xy and monodisperse initial data vq = S(x— 1) [4, 26]: 

iy t = ^2n k (t)S(x-k), n k (t) = k . (8.1) 

k=l ' 6 
A beautiful probabilistic interpretation of this solution in terms of a Poisson- 
Galton- Watson branching process may be found in [4] . The solution is valid 
only for < t < 1. When t = 1, n k (t) only has algebraic decay, and the 
second moment mzit) = oo. Moreover, mass can no longer be conserved for 
t > 1. At a microscopic level, this is commonly ascribed to the formation of 
a cluster of infinite mass (the gel). 

The formal scaling limit of (8.1) is obtained by considering the large k 
limit as t — > 1. By Stirling's approximation k\ ~ V2nke~ k k k , as t —> 1 we 
find 



n k (t) ~ -L= fc -5/2 e *(i-t+lQgt) 



^AT 5 / 2 exp (-fc((l - tf/2 + (1 - tf/3 + ...)) 
2tt 



Let x = k(l — t) 2 and consider the limit k — > oo,t — > 1 such that x is held 
fixed. Thus, we find 

lim (1 - t)' 5 n k (t) = -^Lar^V*/ 2 . 

t— >l,k— >oo V27T 



This shows convergence of the discrete solution to the scaling solution [4] 
n(t,x) = —Lz-^-Ci-t) 2 ^ x€ (0,oo), t€ (-oo.l). (8.2) 

V27T 



We will see below that this scaling solution emerges coherently from the 
scaling solutions to the additive kernel, and is just one of a one-parameter 
family of scaling solutions. 
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8.2 Scaling solutions and weak convergence 

The scaling solutions for the multiplicative kernel can be obtained by our 
knowledge of the scaling solutions to the additive kernel, via a general rela- 
tion between solutions for the two kernels. Recall from Section 2.5 that the 
initial data are normalized so the initial second moment m2(0) = 1. Then 
iri2(t) = (1 — an d the gelation time T ge \ = 1. The second-moment 

probability distribution function 

px J poo 

V(t,x) = J y 2 v t {dy) / J y 2 v t (dy) (8.3) 

is the analogue of M(t, x) for the additive kernel. From (2.22) we see that 
the Laplace transform of V(t, x) is (l — t)d s ip(t, s). We differentiate equation 
(2.24) with respect to s to obtain 

(1 - t)d a il>(t, s) = dM~ log(l -t),s) = u(t, a), (8.4) 

where r(t) := log(l — i) -1 . By consequence, 

V(t,x) = M(r,x), (8.5) 

where M(r, x) is the mass distribution function for the corresponding solu- 
tion with additive kernel. For solutions with densities n(t, x) and n(r, x) for 
the multiplicative and additive kernels respectively, this means 

x 2 n(t, x) = (1 - t)' l xh(T, x). (8.6) 

From this relation we obtain the scaling solutions for the multiplicative 
kernel as described in the introduction. Explicitly, 

n(t, x) = (1 - t)- 1+ Wn„(x(l - t) 1//3 ), (8.7) 

where (3 = p/(l + p) and 

1 ~ (_l)k-l k/3-3 
k=l 

Notice that these scaling solutions do not preserve mass — in fact, all of 
them have infinite mass! Instead, they have a finite second moment for t < 1, 
which blows up as t — ► 1. For < p < 1 the third moment is infinite. When 
p = 1 the scaling solution reduces to the exponentially decaying solution in 
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(8.2) after a trivial scaling. Finally, we note that though we have assumed 
t 6 [0, 1), these solutions are well-defined for t 6 (— oo, 1). 

Theorem 7.1 characterizes the convergence of M(t,Xx) and it is easy 
to adapt to characterize convergence to self-similar form approaching the 
gelation time. 

Theorem 8.1. 1. Suppose there is a rescaling function X(t) — > oo as 
t — ► 1 and a nontrivial probability distribution function V*(x) suc/t 

lim V(t, \(t)x) = V*(x), (8.9) 

at a// points of continuity of V*. T/ien t/iere exists p £ (0, 1] and a 
function L slowly varying at infinity such that 

X 

y 3 u (dy) ~ x 1_p L(x) as x -► oo. (8.10) 

2. Conversely, assume that there exists p £ (0, 1] and a function L slowly 
varying at infinity such that (8.9) holds. Then there is a strictly in- 
creasing rescaling X(t) — > oo suc/t i/iai 

limy(t, A(t)x) = V p (x), 0<x<oo, 

where V p is the second moment distribution function for a scaling so- 
lution, given by V p = M p from (6.4). 

It is worth pointing out explicitly that the domain of attraction of the 
scaling solution in (8.2) includes all initial data with finite second and third 
moments, as well as data whose third moment diverges sufficiently weakly 
(the case p = 1 above). Each of the infinite- mass self-similar solutions, 
however, attracts finite-mass solutions whose third moment diverges at the 
appropriate rate detailed in the theorem. 

The behavior of the rescaling function A(i) and the characterization of 
possibly defective limits can be easily deduced from the corresponding results 
for the additive case that appear in the remarks following Theorem 7.1. 
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